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ABSTRACT 



Context. Gamma-ray bursts (GRBs) were the most energetic events after the Big Bang and they have been observed up to very high 

redshift. Measurements of the chemical abundances are now available for the galaxies hosting such events, which are assumed to 

originate from the explosion of very powerful supernovae ( Type Ib/c), and provide the opportunity to study the nature of these host 

galaxies. 

Aims. We identify the hosts of long GRBs (LGRBs) observed at both low and high redshift to determine whether the hosts are galaxies 

of the same type at different cosmic epochs. 

Methods. We adopt detailed chemical evolution models for galaxies of different morphological type (ellipticals, spirals, irregulars) 

that follow the time evolution of the abundances of several chemical elements (H, He, a-elements, Fe), and compare the results with 

the observed abundances and abundance ratios in galaxies hosting LGRBs. 

Results. We find that the abundances and abundance ratios predicted by models devised for typical irregular galaxies can reproduce 

the abundances of the hosts at both high and low redshift. We also find that the predicted Type Ib/c supernova rate for irregulars is in 

good agreement with observations. Models for spirals and particularly ellipticals fit neither the high-redshift hosts of LGRBs (DLA 

systems) nor the low redshift hosts: in particular, ellipticals cannot possibly be the hosts of gamma-ray bursts at low redshift since 

they exhibit little star formation, and hence no supernovae Ib/c. 

Conclusions. We conclude that the observed abundance and abundance ratios in LGRBs hosts suggest that these hosts are irregular 

galaxies at both high and low redshift, thus demonstra that the host galaxies are the same type of galaxies observed at different ages. 

Key words. Gamma rays: galaxies; ISM. Galaxies: high-redshift. Galaxies: abundances; ISM 



1. Introduction 

Gamma-ray bursts (GRBs) were the most energetic events after 
the Big Bang. They can be detected back to the onset of reioniza- 
tion (e.g., Greiner et al. 2009; Salvaterra et al. 2009; Tanvi r et alj 
2009) due to their brightness in the first few hours after the ex- 



plosion (lLamb & Reichartll2000l) . Long GRBs (LGRBs) (T 90 
2) are transient sources followed by long-lasting afterglows 
emitting energy intensely across the full range of the electro- 
magnetic spectrum. They are effective and informative probes 
of many aspec ts of astrophysi cs, such as cos mic star-formatio n 
rate (SFR) ilKistler et al1l2009t). cos mic dus ylLiang & Lill2009l) . 
galaxy properties (Castro Ceron et al. 2008), and interstellar 
medium (ISM) properties (Levesque et al. 2010a). Whether the 
properties of LGRB progenitors depend or not on metallicity 
is a key issue in understanding the observations. Low metallic- 
ity is indicative of a less evolved galaxy, usually of low stel- 
lar mass, high gas-to-stellar mass ratio, and low luminosity. 
LGRB progenitors a re understood to be low-metallicity massive 
core-collapse s tars (W ooslevl 119931: IWooslev & Hegerf |2006a; 
lYoon et al . 2008) observed in association with supernovae (SNe) 
Ib/c, which in turn originate in massive stars and therefore 
reside in active star-forming regions (e.g. Galamaetal. 1998; 
lHiorfhetal.ll200l iPian et al J 120061: IWooslev & Hegerf l2006b). 
However, this low-metallicity progenitor channel is challenged 
by so-called "dark" LGRB (LGRB051022 and LGRB020819) 



assoc iated with a high-meta llicity environment ( Graha m et alJ 
120091: iLevesque et alfcoiObh . 

Different systems are selected by different methods 
from high to low redshift. Two kinds of LGRB asso- 
ciated systems are selected: LGRB host galaxies at low 
redshift and LGRB-associated damped Lyman-alpha sys- 
tems (LGRB-DLAs) at high redshift. LGRB host galax- 
ies have been studied by severa l authors using ob serva- 
tional data (e.g. [C hary et al ] 120071: iM ichalowski et al. 2008; 
Castro Ceron et al. 2008; Levesque et al. 2010a; Savaglio et al. 



2009; Wolf & Podsiadlowski 2007, among others) and simula- 
tions (e.g. lMaoll2010t IChisari et all|2010tlCamp isi et al. 2009; 
ICalura et alJ l2009t lLapi et all 120081: iNuza et all 120071 among 
others). All these results find that at low redshift (z < 2), the 
typical hosts of LGRB are small, star-forming, low-metallicity 
irregular galaxies. 

At high redshift (z > 2), as for QSO-DLA systems 
(log 
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N H , > 20 cm" 2 ) the LGRB-DLAs are observed 
i Jakobsson et alJl200a iProchaska et alJl2007at iFvnbo et all 
2009). The nature of DLAs remains a matter of debate. The 
DLAs in QSO spectra are assumed to form in the ISM of galax- 
ies located in front of the QSO. The high luminosity of the 
QSOs usually makes it difficult detect the DLA galaxies directly. 
The hypothesi s that these DLAs are the progenitors o f present 
disk galaxies (IWolfe et alJll986HNaab & Ostrikerll2006[) is chal- 
lenged by the low metallicity ( ^ 1/1OOZ ) of DLAs and the 
flat age-metallicity r elation of disk stars (the disk is ra pidly en- 
riched to =* 1/3Zq) dAkerman et alJl2004l: |Pettinill2006l) . Using 
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detailed chemical evolution models, Calur a et ail (120031 [2009) 
investigated the nature of DLAs and suggested that the ma- 
jority of them, including the LGRB-DLAs, may be either spi- 
ral disks observed at large galactocentric distances, irregular 
galaxies such as the LMC, or starburst dwarf irregulars ob- 
served at different times after the last burst of star formation. 
On the other hand, massive elliptical galaxies are unlikely to 
be DLA systems even at very high redshift owing to their in- 
tensive quick enrichment, which produces high [a/Fe] ratios, 
that is inconsistent with the data observed for DLAs. However, 
among the identified high redshift DLA host galaxies, one DLA 
host galaxy is a l uminous Lyman-break galaxy (DLA 2206- 
19 A. lM0ller et al.l 12002). By comparing the luminosity func- 
tions of DLAs and LBGs, IWolfd (120051) concluded that there 
is a significant overlap between the DLA and LBG popula- 
tions. |F^nboetal] ([2008 ) compared the metallicity distributions 
of LGRB-DLAs and LBGs. Their results support the hypothesis 
that LGRB-DLAs could arise from a population of LBGs that are 
not heavily obscured. L BGs are understood to be small young 
star- forming ellipticals (Matteucci & Pipino 2002; Pipino et al. 
2010). Additionally, Zwaan et al.l ( J2005[) showed that in the local 
universe the luminosity distribution of galaxies producing DLAs 
is nearly flat from Mg « -20 to Mb ~ -15, which implies that 
the DLA host galaxies should belong to a complex population 
(a single type population could not span such a wide luminos- 
ity range), it is normally unknown in which region of the galaxy 
the line of si ght from the QSO and LGRBs intersect the ISM 
of the galaxv JProchaska et all d2007bl) argued that LGRB-DLAs 
preferentially probes denser, more depleted, higher metallicity 
gas located in the inner few kpc of the ISM compared to QSO- 
DLAs. This idea is supported by the distribution of Hi column 
densities for GRB absorbers produced using a high resolution 
simulation of galaxy formation dPontzen et al.ll2.Q10b . In princi- 
ple, it should be possible for high-redshift LGRB-DLAs to be a 
phase of low-mass star-forming ellipticals, which is what we test 
in this paper. 

In this article, we use the most recently observed abun- 
dances of the LGRB-associated systems, LGRB-DLAs at high 
redshift and LGRB host galaxies at low redshift, to try to deter- 
mine: i) which galaxies are the hosts of LGRBs? ii) whether the 
LGRB-associated systems represent an evolutionary sequence, 
i.e., whether the same objects at low and high redshift are 
just seen at different phases of their evolution? To do this, we 
adopt updated chemical evolution models for galaxies of differ- 
ent morphological ty pe that reprod uce the properties o f galaxies 
dPipino et al.l d20 1 OJ) for ellipticals. lYin et al.l (12010b) for irreg- 
ulars and spirals), and compare our predicted chemical abun- 
dances with the observed ones. 

The paper is organized as follows: in Sect. [2] we describe our 
galaxy formation scenario; in Sect. [3] we summarize the obser- 
vational constraints; our results and discussions are presented in 
Sect. 21 and our summary and conclusions are drawn in Sect. 
Throughout the paper, we adopt a (0.7, 0.3, 0.7) cosmology. 

2. Galaxy formation scenario 

We now summarize our adopted galaxy formation models for 
galaxies of differen t morphological type. We direct th e reader to 
IPipino et al.l d2010l) (ellipticals) and I Yin et all (1201 Obi) (irregular 
and spirals) for equations and related details. Galaxies of dif- 
ferent morphological type have different star formation histories 
(Matteucci 2001). Different star formation histories produce dif- 
ferent abundance ratios, particularly in terms of [a-/Fe] as a func- 
tion of [Fe/H], while other properties may not show such large 



differences, such as metallicity and stellar mass at high and low 
redshift, respectively (see below). 



2.1. Ellipticals 



The elliptical galaxies are assumed to evolve following an in- 
stantaneous mixing of gas but not according to the instanta- 
neous recycling approximation, i.e., we take into account stel- 
lar lifetimes . The m odel we adopt is similar to that described in 
Pipin o et all (|2010) except that we do not consider the dust pro- 
duction and evolution since her e we c ompare our results only 
with extinction-corrected da ta ([Savaglio 2006; Savagli o et al.l 
2009: ILevesque etaL 2010a), assuming that their dereddenings 



are reliable (but see lLi et al.1 12008). The initial conditions for el- 
lipticals allow formation by either collapse of a gas cloud into 
the potential well of a dark matter halo or, more realistically, 
the merging of several gas clouds. In any case, the timescale 
for both processes, t, should be short (< 0.5 Gyr), so that the 
ellipticals form in a short time. The rapid collapse triggers an 
intense and rapid star-formation process, which can be consid- 
ered as a starburst that lasts until a galactic wind, powered by 
the thermal energy injected by stellar winds and SN (la, Ibc, II) 
explosions, occurs. At that time, the thermal energy is equiv- 
alent to the binding energy of gas, and all the residual gas is 
assumed to be lost. Numerical simulations have demonstrated 
that a reliable assumption is that no more than 20% of the total 
blast-w ave energy of supernovae should be used to thermalize 
the gas (Pipino & Matteucci 2004). The galactic wind develops 
outside-in, and after the wind the star formation stops. Therefore, 
in this picture the outer regions of ellipticals cease the star forma- 
tion before the inner regions. After that time, the galaxies evolve 
passively. An important assumption is that the efficiency of star 
formation, namely the SFR per unit mass of gas, increases with 
galactic mass: this scenario can explain the chemical downsiz - 
ing (for this definition of downsizing, see ICowie et al.l ll996). 
namely that more massive ellipticals have larger [a/Fe] ratios 

in their dominant stellar populations than smaller ellipticals, as 

i — — i j v 

shown by Pipino & Matteucci (2004). In principle, to explain the 

increase in the [a/Fe] ratio one could assume a variable initial 
mass function (IMF) becoming fl atter with incre asing galactic 
mass, or a dry-merger scenario dBell et al.l 120061) . but in both 
of these cases other observations are violated. For instance, 
a top-heavy IMF is inconsis tent with the mass -luminosi ty re- 
lation ship dPadovan i & Matteucci 1 1993b iMatteucci & Padovanil 
Il993l) . while the dry-merger scenario can reproduce the down- 
sizing in [Mg/Fe]-mass relation but not the mass-metallicity re- 
lation and viceversa (see Pipino & Matteucci 2008). One can 
certainly propose a selective metal loss (e.g. Fe should be lost 
more easily than Mg in the more massive galaxies) mechanism 
in any scenario to fit this chemical downsizing, but this as- 
sumption has no a physical basis (see Matt eucci et alJfl998l) . In 
this paper, we adopt a Salpeter (1955) one-slope IMF and the 
same stellar yields as in IPipino et al.l (120101) . In particular, for 
massive stars we adopt the yields suggested by iFrancois et al.l 
(200 4]), for the low and intermediate mass stars the yields of 
van den Hoek & Groenewegen (1997), and for SNe la the yields 
ofllwamoto et al. ( 1999) (model W7). Finally, the assumed SFR 
is a ISchmidtl d!959l) law that is linearly dependent on the gas 
density. 
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2.2. Spirals and irregulars 

Spirals and irregulars are galaxies harboring either recent or ac- 
tive star formation activity. At variance with elliptical galaxies, 
a slow and continuous star formation regime is assumed in spi- 
rals and irregulars. In the irregular galaxies, the efficiency of star 
formation is assumed to be lower than in spirals, an assumption 
that has been proven to be correct. 

Following the work of I Yin et all (|2010a.b), we build a one- 
zone model for irregulars, and assume that the galaxy is built 
up by continuous infall of primordial gas until a given mass is 
accumulated. For the spiral disks, we develop the same model 
as for irregulars but with higher star format ion efficiency (SFE) 
and hi gher luminous mass. As suggested by Calura & Matteuccil 
(2006), the main properties of local galaxies of different mor- 
phological type could be reproduced mainly by decreasing the 
star formation efficiency from early to late types. Therefore, we 
assume for spirals and irregulars a continuous SFR with SFE 
lower than the one used for ellipticals (^3-22 Gyr -1 ), e.g a* 1 
Gyr" 1 for spirals and =* 0.1 - 0.05 Gyr for irregulars. The ir- 
regulars are assumed to assemble by infall of gas on a relatively 
short timescale (t ~ 1 Gyr) and the infall mass is assumed to be 
5 x 1O 9 M for all models, whereas the spiral disks are proposed 
to form by means of a slower accreting process (t ~ 7 Gyr) and 
the total infall mass is 1O U M . Galactic winds triggered by SNe 
are also considered in the same way as for ellipticals. According 
to the different SFE in each model, and consequently different 
SN feedback, the galactic wind develops at different times in 
different models, which then end up with different final stellar 
masses. 

Therefore, the main difference between models for galaxies 
of different morphological type in is the efficiency star forma- 
tion, which is higher in more massive, early-type objects, and 
consequently whether a galactic wind occurs and when. 

3. LGRB data 

The high-redshift LGRB-DLA data are inferred from the rest- 
frame ultraviolet (UV) absorption lines. On the other hand, low- 
redshift LGRB host-galaxy data are measured using the rest- 
frame optical emission lines of H 11 regions. The absorption lines 
in the UV and emission lines in the optical band have little 
overlap in terms of wavelength, making the detection of LGRB- 
DLAs and LGRB host galaxies for the same target particularly 
hard. In this section, we briefly summarize the observations used 
in this paper. 

3.1. Low redshift case: LGRB host galaxies 

The properties of GRB host galaxies were studied by 
Savaglio et al. (2009) using the largest sample studied to date 
(40 LGRBs in total 46 GRBs). Owing to the limitation of mea- 
surements, 89% of the hosts are at z < 1.6. In addition, all 14 
LGRBs considered in this paper for which [O/H] data are avail- 
able, are at z < 1. Ten of them are studied with more de tailed 
optical emission-line diagnostics bv lLevesque et al.l (1201 Pal) . 

3.1.1. Morphology 

The typical morphology of LGRB host galaxies remains unclear. 
The typical observed hosts of LGRB are small, star-forming, 
low-metallicity galaxies. But since there is no reason for LGRBs 
to assume any particular galactic morphology, they should be, 
in principle, observed in galaxies of all morphological types. 



However, by comparing the most faint observed star-forming 
LGRB host galaxies and the rapid evolution of elliptical galax- 
ies, we can deduce that any elliptical hosting a low-redshift 
(z < 2) LGRBs should be small (~ 10 10 M o ) and in its very 
early phase of formatio n. This formation scenario is consistent 
with the model of Pipin o et al.l d2010l) . where we suggested that 
the QSO hosts (~ 10 12 M G ) formed at z ~ 8 and the LBGs 
(~ 5 x 10 10 M Q ) formed at z < 4 implying that at higher red- 
shift (z > 1) one could observe more massive LGRB host ellip- 
tical galaxies. If this scenario were correct, we would observe a 
mass-redshift relation for LGRB hosts. 



3.1.2. Star formation rates 

The SFRs are inferred from either the rest-frame UV (A = 
2800A) continuum luminosity, Ho- flux, or other SFR cali- 
brations (see data references for details). The measured SFRs 
of the adopted samp les span from 0.05 to 36.46 M yr~' 
dSavaglio et al.l I2009I) . Given the complex and uncertain SFR 
measurements, here we adopt only the critical star-forming cri- 
terion, namely that SFR> 0, for the LGRB host galaxies. 



3.1.3. Stellar mass 

The stellar mass can be estimated from observed multi-band 
optical- near-infrared (NIR) photometry by performing spectral 
energy distrib ution (SEP) fitting [j Savaglio et al. 2009). Some 
authors (e.g. lCastro Ceron et al.l l2008: Svensson et al. 2010) de- 
rived the stellar mass only using the rest-frame K-band flux den- 
sity assuming that it is a reliable estimator. Alternatively, the stel- 
lar mass can be estimated using an empirical relati on, when other 
parameters are available (Savaglio 2006). The Savaglio et al. 
(2009) results indicate that the stellar mass of LGRB host 
galaxies are small (maximum M* < 1O H M ). Because of the 
faintness of LGRB host galaxies, and the average observed 
stellar ma ss (< M» >= 1O 93 M ) lower than field galaxies 
(Sav aglio et~atll2009l) . we adopt 5 x 1O 9 M as the typical infall 
mass of our fiducial model for irregulars. On the other hand, we 
consider low mass elliptical galaxies and adopt 1O 1O M as the 



typical mass of our fiducial model of low-mass ellipticals. 



3.1.4. Metallicity 

The metallicity (expressed as log(0/H)+12) was inferred from 
extinction-corrected emission-line fluxes using R23, 03N2- 
metal licity relation, and auroral [O ml A4363 diagnostics by 
Savaglio et al. (2009J) and lLevesque et al.1 (1201 Pal) . The majority 
of the measured oxygen abundances are subsolar. However, two 
super-solar host ga laxies are reported (Leves que et al.l l2010bl 
iGraham et al.l 12009). Those super-solar metallicities challenge 
the assumption of a low metallicity threshold for LGRB pro- 
genitors. As we later show, their hosts may be described by our 
spiral/irregular model with SFE -0.1 (see Fig. [2j. In this case, 
the hosts may be the disks of spirals. 



1 The mass fitting is robust (Shaplev et al. 2005), despite the metal- 
licity, dust extinction, and age being degenerate (see Chen et al. 2010, 
who compared six evolutionary population synthesis (EPS) models and 
concluded that different EPS models are necessary to reproduce differ- 
ent stellar populations). 
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3.2. High redshift case: LGRB-DLAs 

The high-redshift GRB-DLA data are measured by the rest- 
frame UV absorption lines detected in the optical afterglow spec- 
tra. Using the rest-frame UV absorption-line column densities, 
one measures the cold gas-phase (T < 1000 K) abundance, which 
is valid if we assume negligible ionization corrections and that 
the dust extinction corrections were appropriate . The observed 
abundances (e.g. Si, S, Fe) in LGRB-DLA are subsolar, span- 
ning from 1/100 to near solar (Savaglio 2006; Prochaskaet al. 
l2007bHSavaglio et al.l l2009). Since Prochaska et al. (2007]3) did 
not take in to accoun t the dust-depletion correction, we adopted 

the data in lSavag lio (2006J). 

In the literature (e.g. Savaglio 2006), by adopting one reli- 
able element abundance tracer and available dust-depletion cor- 
rections, the metallicity is generally given by that particular el- 
ement abundance. Since the elements are produced by different 
stars with different lifetimes, the metallicity estimated by those 
elements cannot represent the true metallicity ( which is nor- 
mally estimated by [O/H]). Measuring LGRB-DLA metallicity 
is a challenge at z < 1 -6, since the Lya absorption line is in the 
UV. lRau et alJ (}2010) presented an extremely metal-poor LGRB- 
DLA system (Z ~ 1/300 Z ). They measured several chemical 
abundances and suggested that the apparent conflict between the 
high [Si/Fe] and low [O/Fe] ratios found in this object is caused 
by the underestimate of the O column density. 



4. Results and discussion 

The nature of DLAs remains unclear, but many authors have 
suggested that their chemical properties and star-formation his- 
tories can be reproduced by models of irregular galaxies that 
slowly evo lve due to a mild SFR (Bradamante et al. 1998; 
iRomano et alJ2006l : lYin et al. 2010b). Typical elliptical galaxies 
are instead very unlikely to host the DLAs observed at high red- 
shift prim arily because ellipticals form rapidly (in les s than 1 .5 
Gyr, see iPipino & Matteuccill2004t iPipino et allboiOh consum- 
ing and then losing any residual gas so quickly that additional 
star formation is inhibited. Therefore, to observe an elliptical 
when it is still forming stars, one should observe it at very high 
redshift, higher than the typical redshift of DLAs. Moreover, 
if the DLAs represent an early phase of galaxy formation, re- 
gardless of their morphological type, we should observe DLAs 
with different properties ( e.g. different [a/Fe] versus [Fe/H] re- 
lations), which reflect the different morphological galaxy evolu- 
tion history, and this is not the case. 

To verify what has been stated above, we computed chemical 
evolution models for galaxies of different morphological type. In 
Table 1, we summarize some model results: in the first column 
we present the model identification, in the second column the 
total infall mass, in the third column the assumed SFE, in the 
fourth column the predicted present-time gas mass, in the fifth 
column the predicted present-time stellar mass, and in the sixth 
column the predicted present-time (13.7 Gyr) SFR expressed in 
M yr~'. Models from 1 to 4 refer to dwarf irregulars and all as- 
sume an infall mass of 5 x 1O 9 M , whereas model 5 refers to 
an infall mass of 1O U M and the disk of a spiral. Model 6 rep- 
resents a small elliptical with infall mass ~ 1O 1O M and model 
7 a massive elliptical with infall mass of 1O 12 M . We note that 
the predicted SFRs for dwarf irregulars a gree with the lowest 
values estimated by Savagli o et alJ (120091) . while the SFR for 
the spiral disk agrees with the highest observed values, and the 
gas masses agree with the values measured for dwarf irregu- 
lars. No star formation occurs in the elliptical model after the 



galactic wind (occurring before 1.5 Gyr since the beginning of 
star formation), which implies that ellipticals cannot be the ob- 
served nearby (z<l) host galaxies of LGRBs, unless they just 
formed. On the other hand, the outermost regions of galactic 
disks could be associated with DLAs simply because their chem- 
ical and gaseous properties resemble those of irregular galaxies, 
as shown in previous works. 



4.1. Metallicity-redshift relation for LGRB-associated 
systems 

To establish which galaxies host the LGRBs at both low and high 
redshift, we compare our model results with the observed abun- 
dances in LGRB-associated systems (see Figs. Q] [2] [3] and|4]i. 
The model with SFE=1.0 simulates the evolution of a small disk 
of either a spiral galaxy or an irregular galaxy undergoing more 
intense star formation than the average, whereas the SFE=0.1- 
0.05 is typical of dwarf irregular galaxies (see Romano et al. 
2000; lLanfranchi & Matteuccill2003l) . From Figs. Q] and 
|U one can see that the value SFE=1.0 repoduces the LGRB- 
associated systems at neither high nor low redshift, whereas the 
model for irregulars galaxies with SFE=0. 1-0.05 can closely fit 
all the data. On the other hand, values of SFE<0.05 do not fit 
the data. This comparison clearly indicates that the associated 
systems of the studied LGRBs are dwarf irregular galaxies in 
different evolutionary phases. In elliptical galaxies, there is no 
star formation after the galactic wind, therefore they can be re- 
jected as hosts of LGRBs at low redshift. Since ellipticals are 
very rapidly enriched, by a fine-tuning of the formation redshift 
parameter, one can always find a model of elliptical capable of 
reproducing any subsolar log(X/H)+12 data at high redshift, but 
this would be an arbitrary assumption. We therefore need abun- 
dance ratios to test this situation (see Figs. [6] and [7] in the next 
section). 

In Fig. [5] we show again the log(X/H)+12 versus redshift 
plot predicted by a model with SFE=1.0 but an infall mass 
of 1O U M that should represent the disk of a spiral galaxy 
such as the Milky Way (model 5). In this model, we also as- 
sumed that the gas accreted more slowly than dwarf irregu- 
lars. Detailed chemical evolution models for the Milky Way 
(e.g. Chiappini et al. 1997; Boissier & Prantzos 1999; Hou et al. 
2000; Chiappini et al. 2001; Yi n et alJ 120091) predict that the 
Milky Way disk formed in several Gyrs and inside-out. Here, 
we present a simplified one-zone model for a galactic disk that 
predicts average abundances, which helps to verify the behav- 
ior of the age-metallicity relation of a more massive galaxy with 
a higher SFE than typical irregulars. Figure [5] shows that this 
model, as happens for the model with infall mass 5 x 1O 9 M 
and SFE typical of a spiral disk (model 1), cannot fit any of the 
LGRB hosts and therefore should be rejected. We note that we 
assumed a common redshift of galaxy formation for all the mod- 
els (zf = 10) irrespective of the galaxy morphological type be- 
cause old stars are in general found in every galaxy. If galaxies 
formed later, clearly they can not be the hosts of high-redshift 
LGRBs. However, irregulars may be the hosts of low-redshift 
LGRBs if they formed at redshift z=l, although in this case a 
slightly lower SFE would be required. 



4.2. Abundance ratios 

The abundance ratios of chemical elements, which have differ- 
ent timescales of production (such as a-elements and Fe), can 
be used as useful criterion to study the star formation history 
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Table 1. Model predictions at 13 Gyr 



Mod infall mass (M ) SFE (Gyr- 1) Gas mass (lO s M e ) Stellar mass (l0 9 M e ) SFRCMpyr 1 )" 

5*10 9 1 0.079 2.65 0.0079 

5*10" 0.1 0.827 1.34 

5*10 9 0.05 7.34 1.12 

5*10 9 0.01 46.25 0.30 

10" 1 31.25 68.5 

10 10 3 0.0 4.09 

10 12 22 0.0 530 



0.0083 
0.037 
0.046 

31 







in galaxies (see Matteucci 2001). In Figs. [6] and [7] we compare 
the [Si/Fe] and [S/Fe] ratios predicted by irregular and elliptical 
models with t he available dat a for LGRB-DLAs. T he data are 
derived from ISavadiol (120061) and iRau et aT] (1201 Oh . Although 
the error bars in the data are very large, the results imply that 
most of the data cannot be well fitted by the predictions of el- 
liptical models. They are instead reproduced by the predictions 
of irregular models with a SFE=0.05 Gyr 1 , thus confirming the 
results shown in Figs[TJ|2][3] and|4] However, more accurate data 
are necessary in the future to confirm this result. 

In the case both of Si and S, because of the large error bars 
in the observational data, the predictions of the elliptical mod- 
els are marginally acceptable. In general, our predictions indi- 
cate that ellipticals, when forming stars ( a passive evolutionary 
phase is of no interest because is without star formation and SNe 
Ib/c are connected to star formation since they originate in mas- 
sive stars), have higher [a/Fe] ratios than spirals and irregulars 
and that their evolution is much shorter. Hence, if one measure 
were to the [a/Fe] ratios of the gas in star-forming ellipticals at 
high redshift, one should find super-solar values. However, el- 
lipticals are predicted to remain in a phase where their gas has 
a metallicity [Fe/H] in common with those of DLAs (between - 
2.0 and -3.0) for a very short time interval (< 0.1 Gyr), as shown 
in Fig. [8] Therefore, it is very difficult to observe these galax- 
ies in this short time interval, another reason for excluding even 
small ellipticals from being DLA systems. If the DLAs did rep- 
resent a particular evolutionary phase of a galaxy, regardless of 
morphological type, we should then observe DLAs with a wide 
range of properties (e.g. [a/Fe]), which reflecting the different 
star formation histories of galaxies of morphological type. 

Since LGRBs are associated with SNe Ib/c, we also com- 
puted these rates for different galaxies, and compared the pre- 
dicted present-time SN Ib/c rate with the observed one, as shown 
in Fig. [9] To compute this rate, we assumed that the progenitors 
of these SNe are either Wolf-Rayet stars with M,„ > 25M Q or 
massive binary systems in the mas s range 12M < M- m < 2OM , 
as described in detail in iBissaldi et al.l d2007h . For the observa- 
tional SN Ib/c rate, w e assumed the value of the SN Ib/c rate per 
unit mass provided by Mannucci et al. ( 2005 ) for irregular galax- 
ies, which is 0.54+°* SNe (lOOyrT'lO'^Mo. We then multi- 
plied this value by the predicted present-time stellar mass in each 
model (Table 1). The results show that our irregular models pre- 
dict SN Ib/c rates in reasonable agreement with observations. In 
addition, the o bserved GRB redshift distribution peaks around 
z=l (Fig. 1 in Savaglio et al. (2009)) in reasonable agreement 
with our predicted peak of SN Ib/c rate (at z= 1 corresponding 
to a galactic age ~ 5Gyr, see Fig. 9). The sharp break in the SN 
Ib/c rate for ellipticals and irregulars with SFE=0.05 occurs be- 
cause at that time the galactic wind develops (SFE=0.01 is too 



low to develop a wind). Gas is lost through the wind and the 
gas surface density decreases sharply (see Sec.|2]for details). As 
a consequence, the SFR drops and the SN Ib/c rate then also 
drops. 



5. Summary 

We have compared the data of LGRB-associated systems, 
LGRB-DLAs at high redshift, and LGRB host galaxies at low 
redshift, with chemical evolution models for galaxies of differ- 
ent morphological type. We have attempted to answer the fol- 
lowing questions: i) which galaxies are the hosts of LGRBs? ii) 
are the LGRB systems part of an evolutionary sequence, in other 
words are they the same objects at low and high redshift as seen 
in different phases of their evolution? 

Our conclusions can be summarized as follows: 

1. If the observed high-redshift LGRB-DLAs and local LGRB 
host galaxies belonged to an evolutionary sequence, they 
should be irregulars with a common galaxy-formation red- 
shift as high as Zf =10, observed at different phases of their 
evolution. We can fit the majority of these systems using 
slowly evolving irregular models with a star formation effi- 
ciency, SFE=0. 1-0.05. The adopted models were prev iously 
tested on local dwarf irregulars (see I Yin et al.112.0 1 OalTbh and 
they may reproduce their properties. 

2. We tested models with SFE=1.0, wh ich resemble the evo- 
lution of the disk o f the Milky Way jC hiappini et al. 1997; 
Francois et al. 2004; Cescutti et al. 2007). We considered the 
average properties of such a disk, and concluded that they 
fit neither high nor low redshift data since they predict too 
high absolute abundances. We cannot exclude, however, that 
they correspond to the outermost regions of spiral disks, 
since their properties are similar to those of irregulars (see 
lMatteuccietalJll997h . 

3. The rapid chemical enrichment at high redshift and subse- 
quent passive evolution, following the occurrence of a galac- 
tic wind several Gyrs ago, of elliptical galaxies, suggests that 
they cannot be LGRB host galaxies at low redshift and that 
they are very unlikely hosts of LGRB-DLAs even at high 
redshift. The high observed [a/Fe] ratios in ellipticals indi- 
cate in particular that most LGRB-DLAs can be neither a 
phase nor a part of ellipticals. 

4. We have confirmed that the star formation history is the main 
driver of galaxy evolution. In particular, the properties of 
galaxies of different morphological type can be reproduced 
by simply changing the efficiency of star formation. By com- 
paring the observed and predicted chemical properties of 
LGRB hosts, we have shown that they are most closely fit- 
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Fig. 1. 12+log(X/H) as a function of redshift. Predicted lines 
are from the irregular galaxy model with SFE=1, as indicated 
in the figure. High redshift (z > 2) data are LGRB-DLA sys- 
tems and low redshift (z < 1) data are LGRB host galaxies. 
Red pluses, green circles, blue crosses, and black squares are 
relative to O, Si, Fe, and S abundances, respectively. Typical er- 
ror in data is s h own by a vertical line. Data are derived from 
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ted by galaxy models with a lower SFE than either spirals or 
ellipticals, which is typical instead of irregular galaxies. 
5. The elliptical models to which we compared the data refer to 
a very massive elliptical of infall mass 10 12 M Q , which is re- 
produced well by assuming a SFE=22 Gyr -1 , and to a small 
elliptical of infall mass 10 10 M©, which is reproduced well by 
assuming a lower SFE=3 Gyr . These assumptions reflect 
the downsizing in star formation that is necessary to repro- 
duce the majority of the chemical and photometric proper- 
ties of ellipticals, (e.g. downsizing in chemi cal enrichment, 
see iPipino & Matteuccill2004 120061 [2008). According to 
our anti-hierarchical galaxy formation scenario, we should 
with future data observe a mass-redshift relation for LGRB 
hosts. 
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